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Abstract 


The equivalence of the discretized equations resulting from both fluctuation 
splitting and finite volume schemes is demonstrated in one dimension. Scalar 
equations are considered for advection, diffusion, and combined advection/diffu- 
sion. Analysis of systems is performed for the Euler and Navier-Stokes equations 
of gas dynamics. Non-uniform mesh-point distributions are included in the 
analyses. 

Nomenclature 

Geometric and independent variables 

i Computational indicie 

L re f Non-dimensionalizing reference length 

t Edge or element length 

n Normal vector 

r Distance vector 

5 Generalized volume 

t Time 

x Physical coordinate, normalized by L tv j 

T Perimeter of control volume 

ft Generalized integration volume 

Dependent variables 

A Flux Jacobian in conservative variables 

A Flux Jacobian in auxiliary variables 

a Sound speed 

c p Specific heat at constant pressure 

c v Specific heat at constant volume 

E Total energy 

e Internal energy 

F Flux function 

f Numerical flux 

H Total enthalpy 

h Specific enthalpy 

M Subsonic or supersonic matrix dissipation 

P Pressure 

q Heat flow 

R Area-weighted residual 

s Entropy 

T Temperature 

U Conserved variables 

u, v,w Cartesian velocities 
V Primitive variables 
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W Auxiliary variables 

X Matrix of right eigenvectors 

X Right eigenvectors in auxiliary variables 
Z Parameter vector 

k Thermal conductivity 

A Wavespeed 

A Eigenvalue matrix 

p Density 

4> Artificial dissipation function 

(j) Elemental fluctuation 

(j) Fluctuation of auxiliary variables formulation 

(j) ( Elemental artificial dissipation 

<j) f Artificial dissipation in auxiliary variables formulation 

p Coefficient of viscosity 

r Stress component 

Auxiliary symbols 

I Identity matrix 

My Symmetric averaging function 

p, q Arguments of limiter 

P r Prandtl number, P r (air) = 0.72 

5ft Gas constant, 5R(air) — 287 J/(kg-K) 

ft Limiter bound 

f Eigenvalue limiting parameter 

7 Ratio of specific heats, 7 (diatomic) = 1.4 

v Finite element shape' function 

%l) Limiter function 

Operators 

V Gradient 

V 2 Laplacian, V 2 = V ■ V 

A Forward difference, A,-x — Xj+i — 

V Backward difference, = x\ — x*_i 

S Central difference, S(X = ^(xi+i — x,:-i) 

S 2 Second central difference, 6 2 x = AV,x = VA*:r = X{+\ — 2 Xj + x f ;_i 

Acronyms 

COE Contributions from other elements 
LHS Left-hand side 
RHS Right-hand side 

Subscripts 
E Element 
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L Left-hand state 

R Right-hand state 

U Upwind 

2U Second-order upwind 

Superscripts 

i Inviscid 

v Viscous 

Overbars are used to represent cell-average values. Vector symbols indicate 
vectors spanning multiple spatial dimensions. Bold face is used for vectors and 
tensors of systems. Subscripts of variables is short-hand for differentiation. Hats 
denote unit vectors. Tildes denote Roe- averaged quantities. 


Introduction 

Finite volume flux-difference-split, schemes, in particular the Roe scheme[l, 2] 
with a MUSCL[3] second-order extension, are well established for the solution 
of one-dimensional gas dynamics, with textbooks written on the subject [4]. The 
discretization of these schemes on general unstructured domains is well covered 
hv Barth [5]. 

Fluctuation splitting concepts have been introduced for the solution of scalar 
advection problems in two dimensions [6, 7, 8], and are aligned with finite ele- 
ment concepts, as opposed to finite volumes. Notable work has been done by 
Sidilkover[9, 10, 11] to extend fluctuation splitting to the Euler system of equa- 
tions for gas dynamics. 

The current paper systematically establishes the equivalence of the dis- 
cretized equations resulting from both fluctuation splitting and finite volume 
treatments of the Navier-Stokes equations in one dimension on non-uniform 
meshes. The fluctuation splitting development is performed as quadrature over 
discrete elements, in contrast to the usual treatment which resorts to a flux for- 
mula via the divergence theorem. Scalar equations are considered first, covering 
advection, diffusion, and combined advection /diffusion. Then the Euler system 
is considered, equating Roe’s wave-decomposition procedure with Sidilkover’s 
modified auxiliary equations approach. Discretization of viscous and conduc- 
tive terms completes the analysis for the Navier-Stokes equations. 

The equivalence of fluctuation splitting and finite volume in one dimension 
serves as a prelude to multidimensional analysis, where the methods differ. Re- 
sults by Sidilkover[12] suggest there may be definite advantages to fluctuation 
splitting over finite volume for the multidimensional Euler equations. How- 
ever, other researchers [13] have not found an advantage in fluctuation splitting, 
though their treatment of the Euler equations differs significantly from that of 
Sidilkover. 
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Linear data 
representation 
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1 
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2 + 1 


x-axis 


(boundary) 


Control volumes 
and centroids 


Nodes 


Figure 1: One-dimensional finite volume domain. 


Domain 


In one dimension the domain considered is a discretization of the x-axis, either 
with uniform or non-uniform spacing between grid nodes, which are indexed 
by i. Dependent variables are stored as discreet nodal values. 

In a finite volume context a median-dual control volume is constructed about 
each node by defining a cell face halfw T ay between adjacent nodes. This conven- 
tion is depicted in Figure 1, with the cell faces referred to as the ± ^ points. In 
the illustrative case of Figure 1, the 2 + 1 node is a boundary point, and the 
corresponding cell extends only from i + ^ to i + 1. The generalized volume of 
the control cell is, 


Si = x i+ , - Xi _i = 


*E*+1 + + Xj — 1 X ] Xj_ i 


— Six 


(i) 


There is a one-to-one correspondence between the nodes and control volumes. 

Notice that for the nodal distribution depicted, which has non-uniform spac- 
ing, the cell centroids, denoted by x in Figure 1, do not coincide with the nodes. 
Barth [5] states that nodal storage in this case, referred to as mass lumping in a 
finite element context, only alters the time accuracy of finite volume schemes, 
and not the steady state solutions. 

Notice also the piecewise-linear representation of data in the finite volume 
context. Discontinuous jumps in the dependent data are allowed at cell faces. 

Figure 2 depicts the discretization of the domain for the fluctuation splitting 
approach. The data is now continuous and piecewise linear over elements defined 
by the nodes. The centroids of the fluctuation splitting elements are at the 
same locations as the faces of the finite volume cells. No special definition of a 
boundary cell is required. The length of an element is, 


t iy i+ 1 = x ? ;+ 1 -Xi- A iX 


( 2 ) 


4 



Linear data 
representation 


Domain elements 
and centroids 


i - 1 


i + h 


— •— 
* + 1 


x-axis 


(boundary) 


Nodes 


Figure 2: One-dimensional fluctuation splitting domain. 


There are one fewer elements than nodes, and each element is associated with 
two nodes. Correspondingly, each interior node is associated with two elements, 
while each boundary node is associated with only one element. 


Scalar Advection 

A hyperbolic conservation law takes the form, 

U, + V F = 0 


(3) 


where U is the vector of conserved variables and F is the flux of these variables. 
Following Godunov [14], Eqn. 3 can be evaluated in an integral sense, 

[ Utdil = - j V F dft 

Ju Jq 

If the control volume, Q, is fixed in time, then, 

U, dQ = S Q \Jt 


L 


(4) 


(5) 


where the overbar indicates some cell-average value. Using the divergence the- 
orem the flux term can be evaluated as, 


[ V Frffi =<ff-ndT 
Jq Jr 


( 6 ) 


where n is the outward unit normal to the control volume boundary, T. 

The one-dimensional scalar advection problem is obtained from Eqn. 3 as, 


tif -f- F x — 0 


(7) 
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which can be written for a control cell as, 


Sfillf — / F x d£l — Fieft face F right face (8) 

Jn 

Linear Advection 

Linear advection is obtained from Eqn. 7 by choosing F = A u. The advection 
speed, A, is taken to be constant. 

Finite Volume 

Equation 8 is expressed for the finite volume about node i with mass lumping 
to the node as, 


SiUi t — Fj_ ^ fj+l — /i-i /*+£ ~ Ri (®) 

where the numerical flux, /, is a difference expression approximating the exact 
flux function, F. Choosing, 


/<+* = 


F, + Ft 


t+i 


*+2 2 

results in a second-order central difference scheme, 

Ri = —SiF 


( 10 ) 


( 11 ) 


The first-order upwind CIR[15] scheme is obtained by the choice, 


, A + |A| A - |A| 

Ji+k “ 2 Ui 2 U?_l " ] 

Fi + F i+1 |A| A 

= r — A u 

2 2 


(12) 

(13) 


giving, 

Ri = -S t F + S]u (14) 

The first-order upwind is then seen to be the same as a central distribution plus 
an artificial dissipation term, 


*u = hbfu (15) 

Second-order upwind is constructed following the MUSCL concept of van 
Leer[3], where a linear reconstruction is performed on each finite volume. The 
numerical flux of Eqn. 12 is modified to be, 


- A + |A| A — |A| 
/»+* = — + —J-±UR 


(16) 
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where U[_ is the reconstructed conserved variable on the left side of the cell face 
and wr is the reconstructed variable on the right side of the cell face. Following 
Barth [5], a limited reconstruction is performed on each cell as, 

Mface — tij “{- * V face (1^) 

The gradient is evaluated as a central difference, 

(Vu) j= !| (18) 

The limiter function, i/h is employed to provide monotonicity of the solution, 
based upon positivity arguments. The limiter takes the form, 


where, 


^ 7 ± 1 ^7 /V7 \ -* 

P = ~ 1 Q = (VM)i * r i± i 


The more restrictive of the ± choices is used for i /;. Some popular limiters are 
presented in the appendix. 

The discrete numerical flux (Eqn. 16) expands to, 


/ i+ A= 


A + |A| 


'*+* 2 

Fi + F i+ i 


Ui + ip, 


(i.i+l 

25, 


(5 . ?/ ) + 


A-|A| 


t. (19) 


2 4 


/ .. \ 

l«i+i - Vi + 1 2S +I f,+ 1 " J 


^(F+|A| w )-^-<5 I+1 (F 


leading to, 

R t — —SiF + + ^2u 

where the second-order correction is, 


^2U = 


— 1,7 / V^7 — 1 


4 

fi 7 i + 1 




c *i-i(F + |A|u)- ~^(F-|A| 

*^7-1 


»')) 


V’i+i 


^(F+IAH-I^WF-IAIw)) 

Of *J,+ 1 / 


The residual (Eqn. 20) can be rearranged as. 


( 20 ) 


( 21 ) 


Ri = 




5,-i 


r i f V'i+l _ p C 1 

l t I *7,7+1 ^ 1 1 1 ,7 ^ y -C i 


+ 2iPiF i+1 -f 1 , i+l ^ ±i F + , 


5, 


+ 1 


4>2U 


( 22 ) 



where the artificial dissipation is now, 


$2U 


= $U + 


|A| 


—Zi- 1,/-^ ^i-2 + 7T ~ ^i-l 

Di-\ Si 


+ 


V r i-i , „ *l>i+i\ .. ^ ,, 

l,i ' vi,t+l ^ 1 ^ l,i) ^i+1 


5, 


5j 


, V'i+i . 

^M+ 1 rf ^t+2 

•-U+1 


( 23 ) 


On a uniform grid and without limiting, the second-order residual (Eqn. 22) 
reduces to a low-truncation-error central difference minus fourth-order dissipa- 
tion, 


R, = | (-Fi - 2 + 6 F e _! - 6F i+1 + F i+2 )+^ 


-Ui - 2 + 4w,_i 


6u j -h — Ui-i-'i) 

(24) 


Fluctuation splitting 

In the fluctuation splitting framework Eqn. 8 is evaluated over each domain 
element, without recourse to the divergence theorem. The element fluctuation 
is defined as, 

Sq u t ~ fa = - / F x dQ (25) 

J n 

Assuming piecewise linear data, the fluctuation for the (‘ell bounded by x 7 and 
Xi + 1 is evaluated as, 


f Xi +' An 

0i,z+i ~ ^ (M = fi j-j- 1 A — — — A 7 F 

Jxi 


(26) 


The elemental update, the LHS of Eqn. 25, is formed as, 

— ^,»+i ^ g y ~ 2 — ^ Uit + — 0M+i (27) 

Partitioning the fluctuation into halves and distributing equally to the nodes 
yields the elemental update formula, 


^M+l <&,*+ 1 

w . = 

2 ' 2 


^2,2 + 1 

— 1 = 



(28) 


Assembling all the elemental contributions to the nodal updates, it is clear each 
interior node will receive fluctuation signals from the elements adjacent to the 
left and right. The nodal update is formed as the sum of these fluctuation 
contributions, 


U-IA 


u it + 


-«*. = — = s,u„ = ^ ^ (29) 
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or. 


SiU it 


<$>i- 1,?: + 4>i,i + 1 

2 


(30) 


A popular nomenclature convention for Eqns. 28 and 30 is to describe the ele- 
mental distribution formula as, 


SiUi t <— ^ + — 1- COE, Si+iU{+i t *— — h COE (31) 

where COE indicates a sum of similar contributions from other elements joining 
at that node. 

Expanding the nodal update formula (Eqn. 30), 

-V F - A F 

SiU it = f 2 ? = -6, F (32) 

which is the identical central discretization as for finite volume (Eqn. 11). 

An upwind scheme can be constructed by introducing artificial dissipation 
in order to redistribute the fluctuation, 


<£'e = sign(A)tf> E 

The upwind distribution formula becomes, 

^ + COE = 


Si 'll i t 




. rf f~\ <Pi.i+i(l - sign(A)) 


4>e + 4 C0E = <Pu± |(1 +sign(A)) 


(33) 

+ COE 

+ COE (34) 


2 2 
Using the fluctuation definition (Eqn. 26) the nodal update is obtained as, 


Si [Uif — 


(A + |A|)V,u (A - |A|)A,u 


= SiF + 


(35) 


2 2 ’ 2 
which is identical to the first-order upwind discretization for finite volume 
(Eqn. 14). 

A second-order scheme is easily obtained by adding the exact same finite 
volume correction, /? 2 u (Eqn. 21), to the nodal update formula (Eqn. 35). 


Non-linear Advection 

Non-linear advection is obtained from Eqn. 7 by choosing the flux to be, 




% | ^ 
II 

(36) 

Define the Jacobian of the flux, 





-4 = F„ 

(37) 

so that, 





OF 

dF du ^ 



& 

II 

frl 1 

II 

’■ = — A.u x 

on ox 


Equation 

7 may be rearranged in 

non-conservation form, 



Ut + F 

x = u t + Au x — 0 

(38) 
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Finite volume 

Following Roe[l], the analog to the numerical flux of Eqn. 16 becomes, 


-4 l + |-4|j+i 4 r - Hi+i 
fi+k = ^ ^ *-UR 


F l + F r I-4|,+i 


(Ur - U L ) 


( 39 ) 


where ,4 is the conservative linearization, which in this case is, 



ul + 
2 


(40) 


A first-order upwind scheme is obtained using piecewise-constant data, ul = u;, 
and ur = U|._|_ i . A second-order upwind scheme is constructed using the linear 
reconstruction of Eqn. 17. The first-order residual may be written explicitly as, 


Ri = -SiF + 


I "4 1 i+i 





(41) 


Fluctuation splitting 

The elemental fluctuation is 

(pe — — I F x dft — — [ Au x dn 

J n Jn 

Assuming piecewise-linear data Eqn. 42 becomes, 


(42) 


<Pe = -A i+ i A*u = - A {F (43) 

An upwind scheme is created by introducing the artificial dissipation, 


0E = sign(i i+ i)0 E = -|i| i+ iA,t 


(44) 


The distribution formula remains, 


SiU>i t 


Si+iUi+i t f- 


— - ^ + COE 

— + ^ + COE 


(45) 


The nodal update is, 


+ 0j_l,i 01,1+1 - 0+1+1 

Sim, = 2 + 2 

ViF , l-4|,+ i 

= — -YjU — —AjU 

2 2 2 2 

r r, 4 _ l-4|t+5 . 

= -Si F + —2^ A * u 

This is the identical update formula as for finite volume (Eqn. 41). 


( 46 ) 


10 



Expansion shocks 

The discretization of Roe’s scheme allows for unphysical expansion shocks that 
violate the entropy condition. Harten and Hyman [16] proposed a commonly 
used method for perturbing the wavespeeds such that entropy is satisfied and 
expansion shocks are prevented. The correction is applied to any wavespeed 
that can go to zero at a sonic point and takes the form, 

|A| j+ i -f- max 

where the perturbation scale is, 

e — max [o, (A 


|a| h 


r \ 2 


+ e 


(47) 


i+ l - Aj), (A,; + l - Ay + l)j 


(48) 


Scalar Advection/DifFusion 

The governing equation for scalar advection/diffusion problems in one-dirnen- 
sion is, 

u t + F x = Uiu x ) x (49) 


Heat Equation 

Modeling of the viscous RHS in Eqn. 49 begins with a consideration of the heat 
equation, 

u t = (fiu x ) x (50) 

In the finite volume framework one approach to discretizing the viscous term is 
to construct a viscous flux, so that the nodal update becomes. 


SiUi t — (fAU x }jj_}_ (//Ux)j_JL 


(51) 


where, 


{pUx)i+L = Pi+i 


(Vu)j + (Vu)j + i 


(52) 


with the gradients Vu defined by Eqn. 18. This approach leads to a five-point 
stencil. 

An alternative is to use a finite element discretization, which results in a 
three-point stencil. This approach is adopted both by Barth [5] and Anderson 
and Bonhaus[17] in a finite volume context and by Tomaich[18] in a fluctuation 
splitting context. 

A Galerkin finite element discretization, using mass lumping to the nodes, 
is constructed on the fluctuation splitting domain by integrating with the aid 
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of the finite element linear shape function v (see Bickford §4.2.2 [19] or Bathe 
§7.2[20]), 


SiU U = / Vi{fiu x ) x dll 

(53) 

J Q 


Integrating by parts, 


SiU it = Vi(fiu x )\]t\ - / (vi) x (nu x ) dSl 

Jn 

(54) 


The shape function is the linear tent function, and is equal to zero at Xi±\, 
eliminating the first RHS term of Eqn. 54. The remaining term is integrated 
over each element connecting at node i , 


SjUji 



v x u x j.t di } 


(55) 


The dependent variable and shape function gradients are constant over the 
element, and taking the element -average viscosity coefficient the elemental con- 
tributions are, 


Si U; 


A jU 




+1 


£<+* + COE 


A j ii 

Si+iUj+i, e- — - — - — /q+i + COE 


c*,i+i 


The nodal update is written, 


_ A iu _ Vm 

&i u it ~ o Ah+i a Vi-b 


Zu 1+1 


h- l,i 


(56) 


(57) 


Combined Advection and Diffusion 

The combined effects of advection and diffusion in the governing equation 
(Eqn. 49) are treated by discretizing the advection terms as discussed in the 
Scalar Advection section and adding the discretization of the diffusion terms 
from the Heat Equation subsection. Recall, however, that the upwind advec- 
tion discretization includes artificial dissipation, which can mask the physical 
dissipation. 

The best approach for solving discretized advection/diffusion problems, as 
suggested by Barth [21], is to include the maximum of either the physical dif- 
fusion term, as defined by Eqn. 52 or Eqn. 56, or the artificial dissipation, the 
second term of Eqn. 39 for finite volume or ^ in Eqn. 44 for fluctuation splitting. 


Systems 

A hyperbolic conservation law for systems (Eqn. 3) is written in one dimension 

as, 

U, +F X =0 (58) 
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A decomposition of the flux function is sought such that the system can be 
expressed as a decoupled set of advection/diffusion equations. 


Euler Equations 

The one-dimensional Euler equations [22] for perfect gases, suitable for simu- 
lating non-reacting, low-Knudson-number shock-tube flows, are written as a 
conservation law (Eqn. 58) with. 


U = 



(59) 


{ pu 

pu + P 
puH 


(60) 


The Euler equations have a form similar to the non-linear advection problem. 

The total energy and enthalpy are obtained from the internal energy and 
enthalpy, 


E = e + 


u~ 

T 


H = h + 


The energy and enthalpy are related as, 


h - e + 


P 

P 


And the equation of state is, 


P = pe ( 7 - 1) 


(61) 


Finite volume 

The numerical flux remains as in Eqn. 39, 


F|_ + Fr I A|i+^ 




-(Ur-Ul) 


(62) 


Roe[l, 2] constructs the conservative linearization for the |A| f+ i matrix bv 
introducing the parameter vector, 


Z = VP' 


u 

H 


( 63 ) 


The i -h ^ state is taken to be a linear average of the parameter vector, 



Z|_ -f Zr 
2 


13 



Taking the velocity and total enthalpy from the parameter vector, 

H = h 
z 1 Zi 

and defining the Roe-density, 

P = vWr 

the Jacobian matrix is formed as, 


(64) 

(65) 


|A| = X|A|X _1 (66) 

The eigenvalues are, 

A = diag(n, u + a, u - a) (67) 

The right eigenvectors are, 



The product X 1 ! U r — Ul) is evaluated via the characteristic variables, 

2a 2 dp — 2dP ) 
dP + padu > (69) 

dP — padu J 

The sound speed is, 

«'■’ = y = 7(7 - l)c = (7 - 1 )h = ( 7 - 1 )(H - y ) ( 70 ) 

Also note the grouping pdu can be constructed as, 


1 


X (Ur - U L ) = X -1 dU = -r 


2d 2 


pdu — zidz -2 - z 2 dz! (71) 

As for the scalar case, first-order spatial accuracy is obtained by taking the 
right state to be i + 1 and the left, state at i. Higher-order accuracy is obtained 
using gradient reconstruction (Eqn. 17) applied either to each of the conserved 
variables (Eqn. 59) or each of the primitive variables, which are, 


V = 



(72) 


The nodal update is still formed as in Eqn. 9. The residual remains as 
expressed in Eqn. 41, but for systems rather than scalar quantities. 
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Fluctuation splitting 

The Euler flux (Eqn. 60) can be written in terms of the parameter vector, 


F = 


Z X Z 2 


Z2Z3 


Further, the derivative of the flux is, 


dF = 


Z 2 

1=1 z 

1 

0 


Zi 

^z - 2 

Z., 


- Z v 


Z 2 


d Z 


( 73 ) 


(74) 


By assuming a linear variation of the parameter vector on each element, the 
fluctuation is obtained from Eqn. 42 as, 

4>e = - f F x dn = - f FzZ x dn = -Fz&,Z (75) 

Jn Jn 

Deconinck et al [23] show, 

FzAjZ = AA ; U = A,F (76) 

when the Roe-averaged forms (Eqns. 64 and 65) are used to obtain A. 

An upwind scheme is constructed by adding the artificial dissipation, 

4 = -|A| 1+ iA,U (77) 


where | A| is defined in Eqn. 66. Employing the same distribution formula as for 
the scalar advection (Eqn. 45) leads to an update formula analogous to Eqn. 46, 
showing the equivalence between finite volume and fluctuation splitting for the 
one-dimensional Euler equations. 

Before ending the fluctuation splitting discussion, it is desired to frame the 
artificial dissipation in the form, 

<J>E = sign(A I+ i)4>E (78) 

The difficulty lies in defining the matrix sign(A). One expression equates 
Eqns. 66, 75, 76, 77, and 78 to form, 

sign(A)A = |A| = X|A|X* ] 

sign(A) - X|A|X _l A -1 = XIAIA^X" 1 (79) 

Sidilkover[9] offers an alternative to brute force matrix multiplications for 
evaluating Eqn. 79. Introducing the auxiliary variables, W, defined by the 
transformation, 

f/U - UwrfW (80) 
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where, 


with the entropy defined as, 


( ds 
dW = < pdu 


ds = dp 


The Jacobian of the transformation is, 

" 1 0 

U w = u 1 


and its inverse is, 


Un' = 


— u — — h 

2 U 7 -I ^ 2a- 


l-(7" D£r (7-1)^ -(7-1)^ 

—u 1 0 

(7-1) IT -(7-l)« 7-1 


The element fluctuation (Eqn. 75) can be reworked, 

4 > E = -AA f U = -UwU^AUw'U^AiU = -Uu^A.W = V w fe (85) 

where 0 e is the fluctuation as computed for the auxiliary variables, 

0 E - -A&iW (86) 

The flux Jacobian of the auxiliary variable formulation is obtained from the 
conserved flux Jacobian via the similarity transformation, 

A = U^AUw = U^XAX^Uw' = XkX~ l (87) 

so the eigenvalue matrix, A (Eqn. 67), remains unaltered. The right eigenvectors 
are obtained from Eqns. 68 and 84, 

X = U^. 1 X 


* (1) = | o | x (2) = 

The inverse is easily computed to be, 


x {z) = 


1 0 


2^ 2^ 
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The flux Jacobian is evaluated from Eqn. 87, 


A = 


u 0 0 

Owl 
0 a 2 u 


( 90 ) 


which corresponds to the following non-conservative form of the Euler equations, 


St + us x = 0 
put + ufm x + P x =0 

P f + a 2 pu x + uP x =0 (91) 

Having developed an alternative method for obtaining the elemental fluctu- 
ation (Eqn. 85), the artificial dissipation can be addressed (Eqn. 78). 

(j>£ = sign(A)0E = UrvU^?sign(A)Uu — Uu4e (92) 


where, 

0E = U j^. 1 sign ( A ) U w 0 E = sign(^)^ E (93) 

and with the aid of Eqns. 79 and 87, 

sign (A) = U l - v 'sign(A)U, v = U^XIAIA^X^Uw = A’|A|A- 1 A'- 1 (94) 

Using the eigenvalue and eigenvector definitions (Eqns. 67, 88, and 89) sign(A) 
is evaluated to be, 


sign (.4)= 


sign(fx) 0 

0 | (sign(w + d) + sign(u-d)] 

0 f [sign(w + d) - sign(w-d)] 


0 

^ [sign(u + «) - sign(w-d)] 
\ [sign(fi + a) + sign(vi-d)] 


(93) 

By considering two cases, for subsonic and supersonic conditions, Eqn. 95 takes 
on simple forms, 


sign (A) = 


M 8Up if | vx | > a 

M su6 if \u\ < a 


(96) 


where, 


M * up = sign ( vi)/ 


and, 


M suh = 


sign ( w) 0 
0 0 

0 a 


0 

d 

o 


(97) 


(98) 


Navier-Stokes Equations 

The Navier-Stokes equations[24, 25] for the flow of a perfect gas are written in 
one-dimensional conservation law form (Eqn. 58) with U defined in Eqn. 59 and 
the flux defined as, 

F = F‘ - F” (99) 
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where the inviscid flux, F', is the same as the Euler flux (Eqn. 60). The viscous 
flux is. 


F 1 ’ = 


0 

ut xx - q x 


( 100 ) 


Using Stokes’ hypothesis the stress is, 


4 

Fourier’s law for heat flow gives, 


( 101 ) 


q x = -kT x 


( 102 ) 


The thermal conductivity is related to the viscosity through the Prandtl number, 


Pr = 


fiCp 

K 


(103) 


where for air P r = 0.72 [26]. The 
equation of state, 


temperature is obtained from the perfect gas 



(104) 


The inviscid flux is discretized as described in the Euler Equations subsec- 
tion. The contributions from the viscous flux to the nodal update is obtained in 
a Galerkin sense using the system analog to Eqn. 55. No viscous contribution 
is made to the continuity equation. 

Using the linear variation of the parameter vector over an element, the ve- 
locity gradient is locally defined on an element [23, 27], 


(u*)e 




Z2 

T l - 


— (A^2 — uAiZ\) — -r~AiU = ^A iU (105) 

%i z \ p 


where, 


x _2 / y/Pi + y/Pi+l \ ‘ 

P = 2l = ( 2 ) 


(106) 


is called the consistent density average. The viscous contribution to the mo- 
mentum equation can now be expressed, 


L 


4 4 p 

-v I nu x dil = -v x ^i-AiU 


(107) 


The first term of the viscous energy flux is evaluated in a similar manner, 


L 


4 4 _ _ p 

-v x puu x dil = -V x fiu-A,u 
E J 6 p 


(108) 
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The second term requires some manipulation. Begin by defining the temperature 
gradient over an element, 


CTr) E 




where, 



(109) 


( 110 ) 


The heat flow contribution to the viscous flux is then obtained, 


L 


t ' -j kTj (IV\ — 1 K 

E 7^ 


A jP _ A t ^p 

P ~P 


- V 


/t^ /AjF _ 

Pr 7* V P j> ) 


The elemental contributions from the viscous terms is similar to Eqn. 5G. 


SjU 


it 


1 1 1 




ki 


<- 


i+ 1 




o 

if** 

0 




+1 


ff A ' u 




U + 


/ A./* _ \ 

Pr V ^ l ) 


j +COE 

| +COE(112) 


As discussed for the scalar advection /diffusion equations, when solving the 
Navier-Stokes equations the maximum of the viscous contribution to the nodal 
update and the artificial dissipation from the inviscid flux discretization should 
be utilized. When the physical viscous terms are large enough, no artificial 
dissipation is needed. 


Summary 

The equivalence of the discretized equations using both the fluctuation splitting 
and finite volume approaches has been shown for non-uniform one-dimensional 
domains. Advection and diffusion, both separately and together, have been 
considered for scalar equations. For systems, the equivalence of the Roe flux- 
difference-split finite volume scheme and Sidilkover’s fluctuation splitting scheme 
for Euler equations is shown. Finally, viscous diffusion and conduction terms are 
modeled and included, establishing the equivalence for the discretized Navier- 
Stokes equations. The strong link established in one dimension can serve as a 
prelude to multidimensional analysis of fluctuation splitting and finite volume 
as applied to viscous gas dynamics. 
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Appendix — Limiters 

A limiter is a function designed to limit the ratio of two values, satisfying, 

V;(o) = o, ^(i) = i 

A symmetric limiter is defined by, 

Symmetric limiters can also be expressed in terms of symmetric averaging func- 
tions, Af^, obeying, 

q ^\q) = = = p ^(p) (U5) 

A limiter that can achieve a value greater than unity is termed a compressive 
limiter. 

First Order 

The first order limiter, so called because it is usually employed to limit a scheme 
to first order spatial accuracy, is the trivial limiter, 

^ = A/ V , = 0 (116) 


(113) 

(114) 


Minmod 

The minmod limiter[4] is a non-compressive, symmetric limiter defined as, 


4> 


- ) — max(0, min(l, p/q)) 
9 > 


or, 


0 


Vm “ = \ p/q if 


1 


The associated averaging function . is, 
A fy(p,q) = 


pq < 0 

IpI<M 

bl>l«l 


pq <0 

> f IpI<M 

|p|>kl 


(117) 


(118) 


(119) 


The minmod limiter is the non-compressive limit of a generalized 6 limiter 
of Sweby[28]. Minmod is achieved by 0 = 1. The upper limit on 0 is the 
“superbee” limiter, 0 = 2. 





= max[0, min (0p/q, 


1), min(p/g, /?)] 


( 120 ) 
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pq < 0 

ft\p \ < M 

if bl<kl< ft\p\ ( 121 ) 

\q\<\p\< ft\q\ 
ft\q \ < IpI 

pq < 0 

ft\p\<\q\ 

if \p\<\q\< 0\p\ ( 122 ) 

kl < IpI < ft\q I 

ft\q\<\p\ 

A similar, non-symmetric limiter has been proposed by Chakravarthy[29], 

V> = max[0, min (p/q, ft)) (123) 

/ \ f 0 pq < 0 

</, ( - ) = < p/q if \p\ < 0\q\ (124) 

' q ' \ ft |p| > ft\q\ 

The upper bound on this limiter is 1 < 8 < 2. 


i>i p - 


= < 


M(p,q) = { 


0 

ftp/q 

l 

p/q 

ft 

o 

$p 

9 

P 


van Leer 


The van Leer limiter [30] is a symmetric compressive limiter with an upper bound 


of 2. 


fp\ _ y + Ifl _ pq + \pq\ 
L \qj i + £ q 2 + pq 


(125) 


M v = 


pq + HI 
q + p 


(126) 


van Albada 


The van Albada limiter [31] is a symmetric compressive limiter with an upper 
bound of 1.18. 


V' 





pip ± q) 

p 2 + q 1 ' 


pq > 0, zero otherwise 


(127) 


pqjp + q) 

p- + q 2 ' 


pq > 0 


(128) 
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